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Abstract 

One of the key factors behind the rapid evolution of molecular spintronics is the effi- 
cient realization of spin manipulation of organic molecules with a magnetic center. The 
spin state of such molecules may depend crucially on the interaction with the substrate 
on which they are adsorbed. In this letter, we demonstrate, using ab initio density func- 
tional calculations, that the stabilization of a high spin state of an iron porphyrin (FeP) 
molecule can be achieved via a dedicated chemisorption on magnetic substrates of dif- 
ferent species and orientations. It is shown that the strong covalent interaction with the 
substrate increases Fe-N bond lengths in FeP and hence a switching to a high spin state 
(S=2) from a low spin state (S=l) is achieved. A ferromagnetic exchange interaction is 
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established through a direct exchange between Fe and substrate magnetic atoms as well 
as through an indirect exchange via the N atoms in FeP. The mechanism of exchange in- 
teraction is further analyzed by considering structural models constructed from ab initio 
calculations. Finally, we illustrate the possibility of detecting a change in the molecular 
spin state by x-ray magnetic circular dichroism, Raman spectroscopy and spin-polarized 
scanning tunneling microscopy. 

Manipulation and detection of spin states of organic molecules with a magnetic centre are 
important issues in molecular spintronics. The change in the spin state can be achieved by external 
agents like temperature, light, pressure etc. Recently, it was shown ^ that a transition from S=l to 
S=2 spin state of an iron porphyrin (FeP) molecule is possible by strain engineering of a defected 
graphene sheet. The microscopic mechanism behind this change in the spin state is identified as 
the change in Fe-N bond length due to the interaction with strained graphene. The difference in the 
occupancy of molecular orbitals due to different bond lengths is responsible for the change in the 
magnetic moment of the Fe atom at the center of FeP. In this letter, we demonstrate another route 
to realize this effect, but now by dedicated chemisorption on a magnetic metallic substrate, which 
may offer an easier experimental investigation. 

Metalorganic molecules supported on magnetic substrates have attracted a lot of attention in 
the last few years.*^ The primary focus is to understand and manipulate magnetic exchange cou- 
pling between the 3d metal center of a metalorganic (FeP, FePc, CoPc etc.) molecule and a metallic 
magnetic substrate, such as Co, Ni etc. It has been shown that the exchange coupling can be tuned 
to be either ferromagnetic-^ or antiferromagnetic,^ depending on the chemical environment and 
hence on the mechanism of magnetic exchange interaction between Fe and the magnetic atom of 
the substrate. Sophisticated experimental tools, such as element specific x-ray magnetic circular 
dichroism (XMCD)^ and spin-polarized scanning tunneling microscopy (SPSTM)-^ as well as ac- 
curate density functional theory calculations'^ have played an instrumental role in uncovering the 
microscopic mechanisms behind these exotic phenomena. 

The adsorption of metalorganics, e.g., FeP on a metallic magnetic substrate can take place in 
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two distinct ways. One is physisorption where the molecule has a weak van der Waals interaction 
with the substrate without forming any chemical bond and hence, more or less preserves its gas 
phase properties. However, in the scenario of chemisorption, the molecule-substrate interaction is 
much stronger in the presence of well established chemical bonds. ^ A substantial change in the 
geometry and electronic structure of the molecule is expected in this case due to the change in the 
ligand field exerted on the Fe d-orbitals. To be precise, the microscopic reason is the enhancement 
of Fe-N bond lengths ( > 2.04 A), which are larger than those in gas phase or the physisorbed 
molecule. Moreover, we will show below that irrespective of surface orientation, FeP prefers to 
stabihze itself in such a way that it maximizes the chemical bonding between the nitrogen atoms 
of FeP and the underlying substrate atoms. Hence, we propose a general procedure of stabilizing 
metalorganics on magnetic substrates via dedicated chemisorption to introduce a mechanical strain 
on the molecule which leads to a transition from a low to a high spin state. We will also show how 
one can detect the change in the spin states by Raman spectroscopy and SP-STM experiments. 




Figure 1: (Color online) Geometry of a FeP molecule on a Co(OOl) substrate. The shaded central 
part demonstrates the change in the Fe-N bond-distance along with the change in the spin state. A 
smaller Fe-N bond length (< 2.04 A) leads to a smaller exchange splitting with a S=l spin state (top 
right). Both the majority and minority spin channels are thus partially filled. A larger exchange 
splitting is observed for an increased Fe-N bond length (> 2.04 A), where the majority channel 
gets occupied completely and a high spin state (S=2) is obtained (bottom right). 



We have performed first principles density functional calculations using the VASP codeJ^ 
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Plane wave projector augmented wave basis was used in the Perdew Becke Ernzerhof generalized 
gradient approximation (PBE-GGA) for the exchange correlation potential. The plane wave cut 
off energy used was 400 eV. A 3x3x1 Monkhorst Pack k-point set was used for the integration 
in the Brillouin zone. Atoms were relaxed until the Hellman Feynman forces were minimized 
up to 0.01 eV/A. We have used 7x7 lateral supercells of the magnetic substrate for the (001) 
and (111) orientations and a 6x7 supercell was used for the (110) surfaces. In all cases, a slab 
of three layers and 21 A of vacuum length perpendicular to the substrate were considered. For 
geometry optimizations, the atoms in the lowest layer were kept fixed and all the other atoms 
in the slab as well as in the FeP molecule were allowed to relax. To account for the electronic 
correlations in the narrow d- states of Fe in the FeP molecule, we have used GGA+U approach 
where a Coulomb interaction term is added according to the mean field Hubbard U formalism.!^ 
The values of the Coulomb parameter U and the exchange parameter J were chosen to be 4 eV 
and 1 eV respectively as these values correctly reproduce the electronic structure and magnetism 
of FeP in the gas phase.l^For all cases, we have used a semi-empirical approach— to account for 
van der Waals interaction between FeP and the substrates. 

The generality of the change in molecular spin state has been studied by having different types 
of substrates and their orientations, i.e., (100), (110) and (111). Magnetic substrates such as Co 
and Ni with partially filled d-orbitals are chosen to enhance the probability of chemisorption and 
also to allow for a study of the mechanism of magnetic exchange interactions. Different surface 
orientations are chosen to provide different crystallographic environments to the molecule. Ni sur- 
faces with hexagonal (111), rectangular (1 10) and square (100) facets are considered along with a 
Co (100) surface. The (100) surface orientation of Co and Ni provides the most symmetric sub- 
strate (square) for the molecule whereas the rectangular (1 10) or hexagonal (111) surface provides 
a less symmetric environment. This indeed affects the stabilization as well as the geometry of 
the chemisorbed FeP molecule. To be more specific, on a square substrate, Co/Ni atoms interact 
with all four N atoms of FeP in a similar fashion whereas on a rectangular surface, the chemical 
interaction along two perpendicular directions are different. This is also true for the hexagonal 
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surface. We have also considered a configuration where FeP is rotated by 45° with respect to the 
substrate in order to have a less direct interaction between FeP (predominantly with ligand N) and 
the first substrate layer. Even though we obtain a high spin state for all chemisorbed situations, 
chemical stabilization depends on the rotation of the molecule with respect to the substrate, as will 
be discussed below. It is indeed observed that the maximization of the molecule-substrate direct 
chemical interaction increases the stability of adsorption. In all studied systems, can one observe 
that the interaction between N states of the FeP molecule and the surface orbitals of the substrate 
is maximized in the ground state configuration. 

Table 1 : Calculated exchange energies {E'^^ = E^^^ — E^^) and relative energies w.r.t the ground 
state (in parenthesis) in eV are shown for different substrates and their orientations. Note that the 
positive E'^^ indicates that the FM alignment is stable. The configurations denoted by TOP, TOP-R, 
HOLLOW and HOLLOW-R are defined in the text. For the Ni(l 1 1) surface, the TOP configuration 
is the only relevant one. 



Surface 


TOP 


TOP-R 


HOLLOW 


HOLLOW-R 


Ni(llO) 
Ni(lll) 
Ni (001) 
Co (001) 


0.04 (2.5) 

0.13 
0.07 (2.17) 
0.39 (2.49) 


0.09 (0.0) 

0.2 (3.58) 
0.49 (3.6) 


0.43 (1.06) 

0.16 (0.0) 
0.75 (0.0) 


0.13 (1.88) 

0.10(3.13) 
0.32 (4.38) 



The chemisorbed FeP is found to undergo a mechanical strain originating from the molecule- 
surface interaction, primarily by the N-Ni/N-Co direct chemical bonding, which leads to a stretch- 
ing of the Fe-N bond lengths of FeP. The stretching can be symmetric (e.g., for a (001) surface) or 
asymmetric (e.g., for a (110) surface) along different in-plane crystallographic directions depend- 
ing on the orientation of the surface underneath. Nevertheless, all the stretched geometries yield 
a high spin state (S=2) of FeP by having the individual Fe-N bond length exceeding 2.04 A. This 
structural change in FeP has a direct impact on its electronic properties. The four N atoms bonded 
to the central Fe atom (in a Fe^+ state) provides a square planar crystal field. The p-d hybridization 
between N-2p and Fe-3d orbitals shifts the d^2_y2 level higher up in energy due to strong a-type 
hybridization, whereas other orbitals feel a weak ;r-type hybridization and stay lower in energy. 
Along with this covalent contribution, a point charge contribution shifts the d^2_y2 level up due to a 
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strong Coulomb repulsion. Both these parts of the crystal field (CF) are strong when the Fe-N bond 
length is smaller, e.g. in the gas phase. Hence, the CF dominates over the intra atomic exchange 
and a low spin S=l state is observed with 4 spin-up and 2 spin-down electrons. Now as the Fe-N 
bond length is stretched, the CF strength decreases resulting in the lowering of the dj^2_y2 orbital 
([1]). The electrons now occupy Fe-3d orbitals following Hund's first rule by filling up all majority 
spin levels followed by one level in the minority spin channel. Even though we have observed a 
strong bonding between Fe d^i and a surface atom, as the spin state change is mostly governed by 
the shift of the d^2_y2 orbital, CF effects explain the spin state of the chemisorbed FeP molecule. 

Due to chemisorption, the Fe moment has a strong exchange interaction with the substrate 
Co/Ni moments. The magnetic exchange coupling energies {Eex_^AFM _^fm^^ obtained in the 
ground state structures on each surface are given in[T} FM and AFM denote ferromagnetic and an- 
tiferromagnetic alignments between the Fe and substrate moments respectively. It should be noted 
that both direct and indirect superexchange interactions give rise to a ferromagnetic exchange cou- 
pling in these systems. When both Fe and N atoms are on top of surface atoms, a (Co/Ni)-Fe direct 
exchange as well as (Co/Ni)-N-Fe indirect superexchange interactions are operational. This situ- 
ation is denoted as the TOP configuration, for which spin-polarized DOS, energy-resolved charge 
and magnetization densities are shown in |2j The charge density is plotted (inset of [2]^a)) in an 
energy window (shaded part in[2];a)) of 0.23 eV where the hybridization between Fe-N, N-Co and 
Fe-Co orbitals are dominant. The cross section of the total charge density in the (100) plane Qa)) 
indicates a substantial overlap between the orbitals mentioned above. The magnetization densities 
shown in[2]^b) and|2];c) indicate a similar overlap between in-plane and out-of-plane orbitals carry- 
ing magnetic moments. To be precise, the d^2_^i orbital of Fe hybridizes with the orbitals of N in 
the plane of FeP as seen in the spin-down channel Qb)). The out of plane hybridization between 
Fe-Co and N-Co via d^2 and orbitals are observed in the spin-up channel Qc)). Therefore, the 
charge and magnetization densities show the signatures of a direct exchange as well as an indirect 
superexchange interaction between Fe and substrate Co atoms. 

Besides the TOP configuration mentioned above, FeP can be adsorbed in other ways, e.g., TOP- 
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Figure 2: (Color online) (a) Spin-polarized atom and orbital-resolved DOS are shown for FeP on 
Co (001) in the TOP position. In the inset, a charge density cut in the (001) plane is shown for 
an energy window indicated as the shaded region in the DOS plot. Col and Co2 are the atoms 
of the Co substrate, situated just below Fe and N atoms of FeP respectively, (b) Top view of the 
magnetization density isosurface shown for the same energy interval to show the in-plane orbitals. 
(c) Side view of the magnetization density to show the out-of-plane orbitals taking part in the 
exchange interaction between FeP and the Co substrate. 
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R, HOLLOW and HOLLOW-R configurations. For the nomenclatures, the readers are referred 
to the Supplementary Information. If FeP is rotated by 45° angle with respect to the substrate 
(denoted as TOP-R/ HOLLOW-R), N-atoms are no longer placed on top of substrate atoms. In this 
situation the effect of indirect superexchange is minimized and the direct exchange between the Fe 
atom in FeP and the substrate atoms dominates. But the absence of N-substrate bonds pushes the 
molecule further apart from the substrate which in turn, has an effect on the direct exchange. On the 
Co(lOO) surface, FM exchange is favored over AFM by 0.39 eV for the TOP configuration and by 
0.49 eV for TOP-R configuration. For the HOLLOW configuration, the exchange coupling energy 
is 0.69 eV favoring a FM alignment. Therefore, irrespective of chemical species and orientation, 
all surfaces show FM coupling with FeP. Apart from the exchange energies, the relative stabilities 
are also indicated in[T]through the total energies w.r.t. the ground state energy for each surface type 
and orientation. It is clearly observed that the most stable configuration is the HOLLOW position 
for both Co and Ni (001) surfaces whereas the TOP-R configuration is the most stable for a Ni(l 10) 
surface. As mentioned before, stability is primarily dictated by the chemical bonds between the 
atoms in the molecule and in the first surface layer. Also, a difference in behavior in the relative 
stabilities between Co and Ni surface is observed due to the difference in number of J-electrons 
and hence the strength of bonding. 

In order to have a better understanding of the exchange mechanism, we have performed several 
total energy calculations with the structural models extracted from the optimized geometries in our 
full DFT calculations. We have chosen the relevant parts of the optimized geometries to build the 
structural models, which are shown in [3] for the TOP configuration. As we have already stated, 
there are two kinds of possible surface to molecule exchange interactions, viz., (Co/Ni)-Fe direct 
exchange and (Co/Ni)-N-Fe indirect superexchange. [3]^ a) represents the indirect superexchange 
between Fe and Co via N whereas [3]^b) represents the Fe-Co direct exchange. Here Col(Co2) is the 
Co atom just below Fe (N). As mentioned earlier, a rotation by 45° increases the distance between 
Fe and Co. [Sj^c) mimics the situation having a direct exchange but with larger Fe-Co distance. [3fd) 
represents a model, that includes both direct and indirect exchange while [Sj^e) represents Fe-Co 
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Figure 3: Structural models extracted from fully optimized DFT calculations of FeP on different 
substrates indicated in the text. The geometries were not optimized further for these models those 
are used to calculate the total energies for FM and AFM alignments between Fe and Co moments. 
The energy differences between FM and AFM alignments are shown for different models consid- 
ered. 
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direct exchange in the absence of N. ^f) is taken from the optimized HOLLOW configuration. 
In this case, the Fe-N-Co angle is not close to 90° but 75.35° due to the geometry of the surface 
underneath. Based on our calculations, we can conclude that both direct and indirect exchange 
interactions are important to consider with varying importance based on particular situations of 
chemisorption. 

Although the spin moment remains similar in all chemisorbed situations, a deformation in the 
spin density is expected to be quite high for a low symmetry structure, leading to a large value of 
spin dipole moment (T^),^ where (T^) is the expectation value of the z-component of spin-dipole 
operator T. Moreover spin densities differ due to different ligand fields exerted on FeP by different 
surfaces and their orientations. The discussion on spin dipole moment is quite relevant for the 
XMCD measurements where the measured effective moment contains both spin and spin dipole 
moments. We have calculated the spin dipole moments by following the method suggested by 
Freeman et al.^ The calculated values of liT^) are 0.82 /i^, 1.08 /ie and 0.79 Hb for FeP on 
Co(OOl), Ni(lll) and Ni(llO) surfaces respectively. It is clear from the signs and values of (T^) 
that the effective spin moment defined as Sgff = 2{sz) +7{Ty), can be highly affected by (T^). This 
change in the value of Sgff should be detectable in XMCD measurements. 

Besides the spin dipole moments measured in XMCD experiments, we propose two additional 
experimental techniques to detect the change in the spin state of FeP, (i) Raman spectroscopy and 
(ii) spin-polarized STM (SPSTM). We shall first discuss our results for Raman intensities. For 
these calculations, we have considered the optimized geometries of FeP in both chemisorbed and 
physisorbed situations in presence of the substrates and have calculated the vibrational normal 
modes of FeP in those optimized geometries only in the gas phase. The Gaussian 09 program^^^ 
was used at the B3LYP level of theory with the 6-31+G(d,p) basis set. In these calculations, 
the geometries were taken from the optimized structures of FeP on Co(OOl) for chemisorbed and 
physisorbed situations. The calculated wave numbers above 1000 cm^^ were scaled using the 
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Figure 4: (Color online) Calculated Raman intensities for FeP in chemisorbed and physisorbed 
geometries. In the upper panel, Raman intensities are shown for shifts in the range 1 10-450 cm^^ 
whereas the lower panel contains wave numbers in the range 1380-1650 cm^^. In both panels, 
the vibrational modes corresponding to the peaks are denoted along with the corresponding wave 
numbers. 
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correction factor 0.9614^1. Raman Intensities were calculated using the expressiont^i^l 

/(V,-) = /^'Co 

v,(l-e^) 

where / is a scaling factor (1.0 in this case), Si are the calculated Raman activities, Vq is the wave 
number corresponding to the excitation laser line in cm^^ (we used 457.9 nm as it is in the range 
used to measure the resonance Raman spectrum of porphyrins'^, V/ is the calculated vibrational 
frequency, h, c and kg are Planck's constant, speed of light and Boltzmann's constant respectively, 
while T is the room temperature (293 K). 



\ Ef-0.1 eV U I Ef+0.1 evD 




Figure 5: Calculated spin-polarized STM images are shown for bias voltages +/- 0.1 V for FeP 
chemisorbed on (a) Co(OOl), (b) Ni(llO) and (c) Ni(lll) surfaces. STM images for spin-up and 
spin-down channels are shown separately, indicated by up and down arrows respectively. Ef de- 
notes the Fermi energy. 

The Raman spectra of the triplet (S=l) and quintuplet (S=2) spin states differ substantially 
both in the peak intensity and peak position, as evident from |4j There are two spectral windows 
where these differences are significantly observed. First of all, the most intense peaks of the FeP 
(both in the case of the chemisorbed and physisorbed molecule) appear in the range of 1 10 - 450 
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cm^ , shown in the upper panel of |4| The highest band in this spectral window corresponds to a 
breathing vibration of the entire molecule and it is shifted to lower wave numbers in the case of the 
chemisorbed FeP as compared to the physisorbed one by approximately 35 cm ^ A shift of this 
magnitude should easily be observed in Raman experiments.'^ A shift of even larger magnitude 
occurs in the band at 237 cm^ of the physisorbed FeP. It is assigned to a symmetric stretching 
vibration of the entire molecule and it appears 60 cm^ lower in the chemisorbed situation. The 
third major difference in this spectral window is the inversion of the two lower Raman peaks. 
The less intense peaks correspond to the two spectra representing diagonal symmetric stretching 
vibrations of the entire molecule. In the range of 1380-1650 cm^^ (shown in the lower panel of|4]), 
three peaks occur with relatively high intensity in the physisorbed FeP that have correspondence in 
the chemisorbed molecule at much lower wave numbers. For instance, the physisorbed FeP band 
at 1607 cm^^ corresponds to the same C-C-C asymmetric stretching coupled with C-H bending 
that is assigned to the 1425 cm^ chemisorbed FeP band. The other two peaks are assigned to 
in-plane C-C stretching vibrations and C-N-C bendings. Thus, this spectral window would be 
ideal to experimentally detect the switching between physisorbed and chemisorbed states, as the 
disappearance/appearance of these three bands could be easily monitored. 

The other detection method we would like to propose is the SPSTM technique. We have 
simulated STM images using Tersoff-Hamman^I^ approach for spin-up and spin-down channels 
separately. Also, both filled and empty state images are shown for a bias voltage of 0.1 V. A 
similar method has been used by Brede et al^ to simulate SPSTM images for CoPc adsorbed on 
a Fe substrate. From|5} it is clear that for this value of bias voltage, only the molecular states are 
observed and no signature of surface states is seen. Also, one clearly observes distinct differences 
in the features of the images for spin-up and spin-down channels. For example, in the Co(OOl) 
case, the spin-down channel shows a bright intensity in the region of the Fe center, whereas it is 
absent in the spin-down channel for the filled state images. This is true for Ni(l 10) surface as well. 
All these changes observed in the SP-STM images should be useful for detecting the spin state. 

In summary, we have demonstrated by ab initio density functional calculations that a dedicated 
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chemisorption of an FeP molecule on magnetic substrates of different species and orientations can 
switch the spin state of the molecule from S=l to S=2. This change is brought by an increase in the 
Fe-N bond lengths in FeP, induced by the strong covalent interactions between the chemisorbed 
molecule and the substrate atoms. The ground state geometries correspond to the maximization of 
the chemical bonds between N atoms of FeP and the underlying substrate atoms. Moreover, we 
have analyzed the mechanisms of exchange interaction by simple structural models constructed 
from our DFT calculations and demonstrated that both superexchange and direct exchange are 
operational in mediating a ferromagnetic interaction between Fe and the substrate magnetic atoms. 
Finally, we show that it is possible to detect the change in the spin state of FeP due to chemisorption 
by XMCD, Raman spectroscopy and spin-polarized scanning tunneling microscopy. 
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Supplementary Information 

In the Supplementary Information, we provide detailed structural data obtained from our calcu- 
lations. As stated in the main article, stability, spin state and magnetic interaction of FeP with 
the surfaces depend largely on the surface orientation as well as the orientation of the molecule 
itself with respect to the surface underneath. To give an example, the (001) surfaces of Co and Ni 
provide a symmetric crystallographic coordination to FeP, contrary to (1 10) and (111) surfaces. To 
demonstrate the effects of orientations and resulting bonding situations, we have chosen the case 
of Co(OOl) surface. We will demonstrate four possible configurations on top of a (001) surface, 
viz. (i) TOP, (ii) TOP-R, (iii) HOLLOW and (iv) HOLLOW-R. Also, in the following figures, 
chemical bonds less than 2.5A are shown only to make the main content clear. 

(i) TOP: FeP can be chemisorbed in such a way that the Fe center sits exactly on top of a surface 
Co atom making a direct bond between them. The four N's bonded to Fe also orient themselves 
so that each of them makes direct chemical bonds with one Co atom on the surface. We denote 
this particular geometry as the TOP-configuration ([6]). The Fe-N bond length in the FeP-gas phase 
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molecule is 2.0 A where a Co-Co nearest neighbor distance is 2.5 A. To form a direct Fe-Co bond, 
as well as a N-Co bond, which helps in stability, the molecule has to stretch. Therefore, the Fe-N 
bond length is stretched up to 2.13 A M when adsorbed on the (001) surface. 





1.33A 
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(a) 



(b) 



Figure 6: (Color online) FeP molecule on Co (001), TOP position. The gray and green balls 
represent Co substrate atoms. Two different colors are used to present different layers clearly. C 
atoms are shown in steal, while Fe, N and H atoms are represented with blue, orange and yellow 
colors respectively, (a) Top view, where Fe-N bond lengths are shown, (b) Side view, where 
minimum distances between different layers are shown. 



Moreover, a distortion is observed mostly in the surface layer of Co where Co-Co bond length 
is contracted to 2.44 A. Thus mutual stretching and contraction occur during the formation of the 
chemical bonds. As the Fe-N and Co-Co bonds are not exactly the same, an Fe-N-Co bond angle 
larger than 90° is established, 97.44 to 98.16°, to be precise. Fe-Co direct bonding occurs at 2.32 
A even though the minimum distance of the molecule from the surface layer is 2.0 A, as shown in 
[6]^b). The Co atom, just below Fe is pushed down and this also decreases the minimum distance 
between surface Co-layer and the second one down to 1 .33 A while the minimum distance between 
2'"^ and 3'''^ Co layers remains much higher, 1.56 A (|6]l(b)). 

(ii) TOP-R: Another possible configuration with Fe on top of a Co surface atom is denoted as a 
TOP-R configuration, where FeP is rotated by 45 ° around the vertical Fe-Co axis. In this situation, 
N's are placed at hollow sites. This prevents the formation of Co-N direct bonding and hence, the 
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Fe-N-Co indirect exchange path becomes ineffective. The Fe-Co direct bonding distance (2.29 A) 
also changes slightly compared to that in the TOP configuration and so does the minimum FeP- 
substrate distance, which becomes 2.02 A (|7]^b)). The resulting distortion in the Co surface layer is 
also different from the TOP configuration. The Co-Co distance (2.64 A) right below FeP is rather 
stretched than being contracted. The minimum interlayer distance between the topmost and 2"'^ 
layer of Co is also changed to 1.43 A while the same between 2"^ and 3"^ layers changes to 1.68 
A. 

HOLLOW: The distortion in Co layers is considerably less in the HOLLOW configuration 
depicted in [8} In this case, the geometry is stabilized in such a way that the Fe-center does not sit 
on top of any Co atom but on a hollow site, i.e on top of a 2'"^ layer Co atom. As there is no direct 
bonding between Fe and Co, the surface Co layer is not much affected. The Fe-Co bond length 
is much larger here whereas a direct bond exists between N's and surface Co atoms. In this case, 
next nearest neighbor Co-Co distance (3.60 A) should be compared with N-Fe-N distance (4.25 A) 
(jSj^a)). This forces the Fe-N-Co bond angle to be 75.35 °, which is much smaller than 90 °. The 
minimum distance between FeP and the surface layer is 1.78 A. The minimum distance between 
-^st_2nd 2f^d_yd ^^ygj-g j^qj differ much (1.62 Aand 1.66 A, respectively), as shown injsj^b). 




(a) (b) 

Figure 7: (Color online) TOP-R configuration, (a) Top view, where Fe-N bond lengths are shown, 
(b) Side view, where minimum distances between different layers are shown. 
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Figure 8: (Color online) HOLLOW configuration, (a) Top view, where Fe-N bond lengths are 
shown, (b) Side view, where minimum distances between different layers are shown. 




Figure 9: (Color online) HOLLOW-R configuration, (a) Top view, where Fe-N bond lengths are 
shown, (b) Side view, where minimum distances between layers are shown. 
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HOLLOW-R: If we rotate FeP with respect to the HOLLOW configuration by 45 ° about the 
vertical axis (denoted as HOLLOW-R), Co-N direct bonds will not form and as in the previous 
structure, Fe-Co bond is longer compared to that in TOP/TOP-R configuration. Unlike the previous 
one, the indirect Fe-N-Co exchange path is ineffective here. The Fe-N bond length is 2. 10 A (|9]^a)), 
which is slightly less compared to other structures but is enough to change the spin state of Fe. The 
minimum distance between FeP and surface Co-layer is 1 .89 A. As the distortion is not much in the 
first Co layer, the minimum distances between l^'^'-2"'^ and 2"'^-3"^ layers are pretty much similar, 
1.59 A and 1.66 A respectively Qb)). 
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